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Electric fields in solar magnetic structures due to gradient 
driven instabilities: heating and acceleration of particles 
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ABSTRACT 

The electrostatic instabilities driven by the gradients of the density, temperature and 
magnetic field, are discussed in their application to solar magnetic structures. Strongly 
growing modes are found for some typical plasma parameters. These instabilities i) im- 
ply the presence of electric fields that can accelerate the plasma particles in both 
perpendicular and parallel directions with respect to the magnetic field vector, and 
ii) can stochastically heat ions. The perpendicular acceleration is to the leading or- 
der determined by the E x B-drift acting equally on both ions and electrons, while 
the parallel acceleration is most effective on electrons. The experimentally confirmed 
stochastic heating is shown to act mainly in the direction perpendicular to the mag- 
netic field vector and acts stronger on heavier ions. The energy release rate and heating 
may exceed for several orders of magnitude the value accepted as necessary for a self- 
sustained heating in the solar corona. The energy source for both the acceleration and 
the heating is stored in the mentioned background gradients. 

Key words: Sun: atmosphere - Sun: oscillations. 



1 INTRODUCTION 

Very strong electric fields have been reported in the solar 
atmosphere in many studies in the pa st. The ma gnitude of 
such fields may amoun t to 7T0 4 V/m ([Pavisll 19771 ) and even 
up to 1.3 ■ 10 5 V/m <|Zhang fe Smarttl Il98fj ). It is widely 
believed that these electric fields appear in the process of 
magnetic reconnection. The magnetic reconnection itself im- 
plies a considerable change in the magnetic field topology, 
and it is also used in describing rapid energy releases and 
'energization' of plasma particles, i.e. their heating and ac- 
celeration. However, at least in some cases 



iMavfield fc Chapmanlll98ll : iPudovkin et al.l 



Jansscns 1972; 



1998) such en- 



ergy release events seem to appear without a measurable 
change of the magnetic energy and configuration and, hence, 
a different description and source of such events may be 
needed. In the present work, we show that such an alterna- 
tive source naturally exists, and that it may drive extremely 
strong electrostatic instabilities and produce electric fields of 
the order mentioned above. The energy for such instabilities 
is provided by the omnipresent gradients of the background 
plasma parameters (the density, temperature, and magnetic 
field). However, the physics and the features of this source, 
and the mechanism by which its energy is transferred into 
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these instabilities are beyond the standard magnetohydro- 
dynamics (MHD) model. In order to describe it, a multi- 
component (fluid or kinetic) theory is needed. 

Typical studies of oscillations and instabilities in the so- 
lar magnetic field structures imply a cylindric geometry with 
parameters having different values inside and outside of the 
structure, yet having at the same time some constant values 
in the two separate domains. Moreover, the usually curved 
structures (e.g. magnetic loops) are often flattened out in 
the modeling and the curvature effects are thus neglected 
in the simplified models. In reality however, the plasma pa- 
rameters (the density, temperature and magnetic field mag- 
nitude) may vary both in the axial (parallel to the mag- 
netic field vector) and in the radial (perpendicular) direc- 
tion. The axial variation is typically on a much larger char- 
acteristic spatial scale in comparison to the radial variation, 
and in some cases can be omitted. The presence of these ra- 
dial inhomogeneities in the background of some accidental 
fluctuations implies a source of energy for certain types of 
instabilities. These background (equilibrium) gradients also 
imply the equilibrium drift velocities, and the associated in- 
stabilities are usually terme d as reactive drift in stabilities, 
first predicted long ago by iRudakov fc Sagdeevl l|l96ll ) . A 
drift wave driven by the density gradient is typically grow- 
ing due to the electro n thermal effects in both collisio nal and 
collision-less regimes (|Vranies fc Poedtsl l2006, 2009aI3). The 
former is well described within the two-component fluid the- 
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ory, the latter however, is a strictly kinetic effect. In both 
cases the ions play a stabilizing role, and in some situations 
they may even impose a threshold for the instability. 

However, in the case of hot ions and in the presence of 
both the density and temperature gradients, the above men- 
tioned reactive instability is termed as ^-instability, where 
now the ions play a crucial destabilizing role. Here, r\i = 
L n /L T , and L n = (dno/dx/no) , L T = (dTo/dx/To) -1 are 
the characteristic inhomogeneity scale-lengths of the equi- 
librium quantities that are here, and further in the text, 
denoted by the index 0. The coordinate x in the present 
local analysis is used to describe the changes in the radial 
(perpendicular) direction. The background magnetic field is 
typically also inhomogeneous (i.e., with a curvature and a 
gradient in the perpendicular direction), and this may be 
described by yet another characteristic scale- length L B = 
(dBo I dx I Bq)~ . The interplay of these three gradients de- 
termines the behavior of low frequency (ui <C Q,i = qiBo/rrii, 
qi = Zie is the ion charge) electrostatic oscillations and 
instabilities. This will be demonstrated in the forthcoming 
text. 



2 ELECTROSTATIC INSTABILITY IN AN 
ADVANCED FLUID MODEL 

We apply the two-fluid model deve loped in numerous 
works related to labora tory plasmas ( Nil sson et al.l Il990l ; 
INilsson fc Weilandl ll 995), A systematic presentation of the 
theory that has been successfully used in the past in the 
predicti on of tran s port p rocesses in tokamak plasmas can be 
found in lWeilandl (j200Cj). There, one can also see a complete 
agreement between this advanced fluid model and the kinetic 
theory (that is one of the reasons for the term 'advanced' 
used here) . Part of the basic th eory of the drift wa v e app lied 
to the solar plasmas is given in IVranies fc Poedtsl ibOOrl). It 
has also been used very recently (|Vranies fc Poedtdl2009al lbT) 
in order to explain some essential properties of the coronal 
heating mechanism. The theory is described in detail in the 
references mentioned above. For completeness, we shall pro- 
vide here a general description of the derivations, emphasiz- 
ing some most important features of the model and provid- 
ing explanations for the assumptions used in the procedure. 
The present analysis is restricted to the electrostatic limit. 
Note however, that th e theory works well also in the full 
electromagnetic limit l|Andersson fc Weilandl [i98Sf ). where 
it can be used e.g. for describing the ballooning instabilities. 
This domain also can be of great importance for the solar 
plasma as it may provide a triggering mechanism for abrupt 
changes in the magnetic field topology, e.g. in processes like 
magnetic reconnection and Coronal Mass Ejections (CMEs). 

The presence of hot ions (typical for the solar atmo- 
sphere) and the background temperature gradient (that is 
expected in solar magnetic configurations) implies, first of 
all, the necessity of including their full thermal response (the 
pressure and the gyro- viscosity collision-less stress tensor) in 
the momentum equation, and, second, the use of the ion en- 
ergy equation in the mathematical model. For the present 
purpose, the later comprises the diamagnet ic heat flow term 
only, and can be written as (|Weilandl l200C0 




Figure 1. Contour plot of F(L n , L T , Lb) = Vi ~ Vthi f° r k y p s 
0.1. The positive values correspond to the instability. 
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Here, Ti is in energy units, q„i is the diamagnetic 
heat flux, and e\\ = B/B. The given form of q*i 
can be obtained directly from the drift-kinetic theory. 
In the case of different temperatures ( pressures) in the 
two directions (|Mondt fc Weilandl Il99ll ) it is to be re- 
placed with \pix/(mi£li)]e\\ x V(2T l± + T q /2) + [(p q - 
Pi±)/(fii^i)]Ti||e|| x (en • V)ey. A detailed analysis of the 
temperature anisotro py effect s on th e gradient driven insta- 
bility is performed bv iMondtl (|l996l ). 

The magnetic field is inhomogeneous in the general case. 
This implies that, in the continuity equation, the contribu- 
tion of the diamagnetic drift to the ion flux does not vanish, 
and the appropriate linearized term is 

V • (rnvti) = ^v bi ■ Vp t i / 0. 

For the same reason we have also V ■ «e / 0, where v E is 
the E x B-drift. Hence, an additional magnetic drift veloc- 
ity Vbi appears in the description of the ion motion. We 
use standard notation from the drift wave theory where 
v*i — e 1 1 x X7pi/(qiriiB). In the equations above we have 

V • q»i = -5n t v« t X7Ti/2 + 5rnv bi ■ VTi/2. The first (non- 
curvature) part in this expression cancels out in the proce- 
dure of calculating V • Vi in the ion continuity equation. The 
second term comprises the ion magnetic drift, which in the 
general case is the sum of the curvature and the gradient-B 
drifts 
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i + Vi ■ V^j Ti +piV ■ Vi = -V ■ q,i, 



(1) 



Here, (eii ■ V)ey = —R/R 2 and R denotes the radius of the 
curvature of the magnetic field, while the two velocities are 
in general case different v 2 = Ti/rrii, v± = 2Ti/m,i. In what 
follows, we sha ll use the exp ression for the effective total 
curvature drift (Weiland 2000) Vbi ^ [2Tio/(qiBo)]e\\ x (ey ■ 
V)e||. 

The perturbed ion perpendicular velocity can be ob- 
tained from the ion momentum equation by applying the 
vector product eii x yielding 
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(2) 



Here, v E and v,i are already denned above, the third term 
is the ion polarization drift, and v„i denotes the drift due 
to the ion stress tensor effects. For small accidental per- 
turbations propagating predominantly in the perpendicular 
direction ~ exp[— iujt + ik y y + ik z z], \k y \ S> \k z \, \ ui\ <S fii, 
the linearized ion energy equation yields ((Weiland 20001) 
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Here, the terms ui*j , u>bj are the product of the perpendicular 
wave number component k y and the diamagnetic and mag- 
netic drifts, respectively. Note that ui, e = k y v* e — —to;*.;, 
t = Teo/Tio- Using Eq. ([2]) in the ion continuity, th e ion 
density perturbation can be written as |Weilandll2OO0h : 
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The electron parallel dynamics yields just the Boltzmann 
distribution for the electron number density, and using the 
quasi-neutrality condition, one then obtains the dispersion 
equation 



ff (l + k 2 y p 2 s ) + n 
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Here, Q = w/w» e , Ps = c s /fii, = T e0 /mi, n 6 = uj bi /uj ti = 
L„/L B , and 0Ji, e = —tojm- It is seen that, without the mag- 
netic field inhomogeneity, Eq. yields only the standard 
drift mode driven by the density gradient, with lo ~ 1/L n . 
The magnetic drift terms lou are responsible for the appear- 
ance of the additional plasma mode, while all three gradi- 
ents together a re responsible for the instability. As shown by 
IWeilandl (|2000l ). Eq. © can be solved analytically yielding 
the approximate growth-rate (normalized to uj, e ) 
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It is seen that for rji > n t h there will be an instabil- 
ity. For a fixed r, the stability conditions are completely 
determined by the three characteristic inhomogeneity scale- 
lengths L n ,L T , L B - This is demonstrated in Fig. 1 where 
we give the contour plot of F(L n , L T , L B ) — rji — rj t h for 
k y p s = 0.1 and for r = T e o/Tio = 1- The positive lines 
denote the values for which the gradient-driven instabil- 
ity takes place. In application to the solar atmosphere with 
T e0 = T l0 = 10 6 K, and assuming B = 3 ■ 10" 2 T, for hy- 
drogen ions we have p s — 0.032 m and, hence, the condition 
k y p s =0.1 implies the perpendicular wavelength \ y = 2 m. 
The assumption of the nearly perpendicular perturbed ion 
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Figure 2. Two real solutions (normalized to u)* e ) in the case 
rji = 1 and k y p s =0.1. 



motion for such a short perpendicular wave-length in fact 
implies a flute-like mode that is very elongated along the 
magnetic field vector, with a parallel wave-length that is 
measured in hundreds of kilometers. Such strongly elon- 
gated modes (i.e. k± 3> fcii) are easily excited under lab- 
oratory conditions (e.g. in a tokamak plasma) in spite of the 
rather limited scales in the parallel direction. In the solar 
magnetic structures with naturally drastic differences in the 
perpendicular and parallel scale-lengths, their excitation is 
expected to be even more efficient. 

The second order dispersion equation Q is solved nu- 
merically and some results are presented in Figs. 2-5. In 
accordance with Fig. 1, in Fig. 2 for rji = 1 we have two 
real solutions for the wave frequency, one essentially due to 
the density gradient and the other due to the magnetic field 
gradient. The initially positive solution changes the sign for 
rjt — 0.4 and then both solutions propagate in the direction 
of the ion diamagnetic drift. Note that for the parameters 
used above and for k y p s — 0.1, the normalized frequency 
Q ~ 1 Hz here implies L n = L T ~ 10 3 m. 

In Fig. 3, the case rji = 3 is presented. Here, both ini- 
tial solutions are positive and real for small values of rji,- 
They merge at around nt, ~ 0.09, yielding a pair of complex- 
conjugate solutions with the real part becoming negative 
for rjt > 0.25. The instability vanishes for nt > 2.04 when 
two negative real solutions appear. The mode is particularly 
strongly growing (|-y| > |fl r [) in the range nj £ (0.15, 0.7). 

Considering rji — 5, in Fig. 4 we present the same mode 
behavior as in Fig. 3. Such a larger rji value yields a growth- 
rate that is larger by about a factor 2 and the instability 
range in terms of rjb is widened to rjb G (0.05, 3.7). 

In Fig. 5, for a fixed rjb = 1 and k y p s — 0.1, the fre- 
quency is calculated in terms of rji. The two negative real 
solutions (c.f. Fig. 2) merge for rji ^ 1.78 yielding a pair of 
complex conjugate solutions with a very weakly decreasing 
real part (£l r — —1.67 at rji = 1.78 and f2 r = —1.71 at 
th = 10). 

Note that in all these cases the used values of the ratios 
rji = L n /L T , rjb = L n /L B in fact imply a wide range of 
possible values for L n , L T ,L B - Since u) r = fi r u* e ~ l/L n , 
this also implies a wide range of possible frequencies. 
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Figure 3. The solutions of Eq. {2) for rjj = 3 showing two real 
solutions for small r/i, and for % > 2.04. The dashed line is the 
growth-rate for the complex-conjugate solutions in between. 




Figure 4. The solutions of Eq. 1(4} for r\i = 5. 
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Figure 5. The solutions of Eq. Q in terms of for r;^ = 1. 



Figure 6. The threshold value of the potential <j>, Eq. , in 
terms of the perpendicular wave-length, above which the stochas- 
tic heating takes place. Above the dotted line the parallel electric 
field exceeds the Dreicer value, resulting in a simultaneous accel- 
eration of particles. 



3 ELECTRIC FIELD 

3.1 Acceleration of plasma particles 

The electrostatic instability discussed here implies an elec- 
tric field varying in time and space, and having very dif- 
ferent components in the parallel and perpendicular direc- 
tions. Such an electric field can accelerate plasma particles 
in both perpendicular and parallel directions with respect to 
the magnetic field vector. An acceleration always exists in 
the presence of electrostatic perturbations, yet typically it 
is sporadic and acts mainly on the small amount of particles 
from the far tail in the distribution function. The critical 
value of such an electric field, above which the bulk elec- 
tron runaway effect takes place, in a fully ionized plasma is 
|Dreicerlll959l ) E d = eL e i/{^e \ d ). Here, L ei = log(A d /6) 
is the Coulomb logarithm, Ad = A de A (ii /(A de + A 2 ^) 1 ' 2 is the 
plasma Debye radius, Adj = v T j/io P j, v T j,u P j are, respec- 
tively, the thermal velocity and the plasma frequency of the 
j species, and b — [e 2 /[127teo (T e + Ti)] is the impact param- 
eter for electron-ion collisions. For the parameters used in 
the previous text we have L e ; = 19, Ad = 0.0005 m, and 
the Dreicer field is 0.11 V/m. Assuming the parallel wave- 
length of about 100(500) km, the amplitude of the elec- 
trostatic potential <j) necessary to achieve the Dreicer value 
is about 1.8(9) KV. However, in the perpendicular direc- 
tion this same potential gives the electric field k y (j> that is 
around E± — 5.7(29) KV/m. For the parallel wave-length of 
about 1000 km we would have <f> = 18 KV and consequently 
E± = 57 KV/m. These estimates are for the number density 
no = 10 16 m" 3 . 

Taking the number density one order of magnitude 
larger no = 10 17 m" 3 yields E d = 1 V/m. The electric field 
corresponding to this value in the parallel direction would 
(for the three given parallel wavelengths) in the perpendicu- 
lar direction have the magnitude of 54, 270, and 540 kV/m, 
respectively. 

The time needed for the perturbations to achieve such 
values can be estimated from the previously calculated 
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Table 1. Stochastic heating for hydrogen (and helium 
in brackets) for several values of the perpendicular 
wavelength and the wave amplitude 0. The assumed 
starting temperature is 10 6 K. 



\ y [m] 


<t> [V] 


a 






0.1 


50 


2.3 (9.1) 


1.5 


10 6 (2.6 ■ 10 6 ) 


0.2 


100 


1.1 (4.6) 


3 


10 6 (3.5 ■ 10 6 ) 


1 


3000 


1.4 (5.5) 


9 


10 7 (1.2 ■ 10 8 ) 


1 


5000 


2.3 (9.1) 


1.5 


10 s (2.6 ■ 10 8 ) 


2 


10000 


1.1 (4.6) 


3.1 


10 8 (3.5 ■ 10 8 ) 



growth-rates. Taking as an example Fig. 3, the maximum 
growth-rate is 7/u, e ~ 1. Taking L n — 10 3 m, Tj = T e = 
T = 10 6 K, B = 3 ■ 10" 2 T, for X y = 2 m we have 
oj* e = 9 Hz. Assuming some small starting value of the elec- 
trostatic potential ty, the growth time till it gets some value 
4>! is t g ~ log(^i/vf)/7. Taking e*/(/tTi) = 0.01 this yields 
$ = 0.86 V. Hence, the value cj>i — 9 KV discussed above is 
achieved within t g ~ 1 s. Observe that for L n = 10 4 m, this 
growth time becomes 10 seconds. 

We conclude that the presented instability can yield 
the extremely la rge values of the electric field repo rted in 
the observations (|Davislll977l ; IZhang fe Smarttlll986D within 
seconds. The electric field generated in such a way may ex- 
ceed the Dreicer value, and consequently an acceleration 
of the bulk plasma may take place. In the direction per- 
pendicular to the magnetic field vector, the particles are 
subject to the leading order E x B-drift v E ~ Ex/ Bo, 
which is the same for both electrons and ions. Note that 
for E± = 20 KV/m and the earlier assumed magnetic field, 
it is of the order of 700 km/s. Hence, because of such short 
perpendicular wave-lengths, these are very short-scale per- 
pendicular plasma fluxes acting on the plasma as a whole. 

The nature of the instability is that in the perpendic- 
ular direction (here along the a;-axis) it is localized in the 
area of the maximum gradients e.g. around some value xq. 
On the other hand, the perpendicular component of the per- 
turbed electric field is in the indirection (corresponding to 
the azimuthal direction in the cylindric geometry), so that 
due to the E x B-drift the plasma fluxes will be in the re- 
direction. Hence, the starting regular density profile in the 
^-direction will be modified and density condensations will 
appear one after another as one moves in the {/-direction, left 
and right of the position xq. In a realistic cylindric geome- 
try and taking the mode behavior in the axial direction into 
account, this would imply the formation of braided density 
structures twisted along the cylinder. However, in the limit 
of large potential amplitudes the initial plasma configuration 
can simply be destroyed. 

In the parallel direct ion the velocity ijBittencourtll 19951 : 
IVranies fc Poedtsl l2009"bl ) is v jz {t) ~ E z /[mj(k z v - u r )]. 
Hence, the acceleration is proportional to mj 1 and acts 
mainly on electrons. It is selective in the sense that par- 
ticularly strong acceleration is experienced by resonant elec- 
trons having the starting velocity vq satisfying the condition 
vq z = u) r /k z . More detai ls on that issue may be found in 
iFletcher fc Hudson! (|2008f l and lVranies fc Poedtsl |2009bl ). 



3.2 Quasi-static purely growing instability 

The instability described in the previous text implies the 
presence of a range of values for the three parameters 
L n , L T , L B for which |fl r | -C 7- This is due to the demon- 
strated change of the mode direction, during which the dis- 
persion lines intersect with the 7/(,-axis so that cu r — 0. In 
Figs. 3 and 4 this is around rjb — 0.25. The corresponding 
growth-rate for the two cases is 7 ~ 0.6 and 0.9, respectively. 

This implies an almost purely growing, quasi-static elec- 
tric field, yet spatially varying and with its amplitude de- 
termined by the the two mode numbers k VtZ . The above de- 
scribed acceleration will remain similar, yet the important 
difference is that the mode is practically non-propagating 
and the spatial variation of the acceleration will become 
much more pronounced and the mentioned braiding more 
effective. 

3.3 Stochastic heating 

The polarization drift, i.e. the third term in Eq. starts to 
play an important role for a large enough wave amplitude. 
In this case, the motion of a particle becomes stochastic 
and consequently he ating takes pl ace. Details of this pro- 
cess can be found in iBellanl (|2006). It turns out that for a 
large enough wave amplitude the standard iterative proce- 
dure, which is behind Eq. ((2j), is not valid any more, and 
the same holds for the particle representation by its gyro- 
center. Instead, one is supposed to describe the actual par- 
ticle motion by writing the particle momentum equation for 
its motion in the wave-field. This has been described in de - 
tail and experimen t ally v erified in McChes nev et all (|l987l ) 
and ISanders et all fl998). It is shown that the stochastic 
heating takes place provided that 

, 2 2 

a = k yPi ■ 7f~ > !- (5) 

Here, pi = v T i/0.i. The condition ([SJ in fact implies that in 
this regime the ion displacement due to the polarization drift 
is comparable to the perpendicular wavelength. This is be- 
cause v p i — (e z x dvj ± /di)/Qj, and Vj± is the leading order 
E x B-drift, so that |McChesnev et al.|[l987l ) v p i ~ auj/k y 
and the perpendicular displacement due to the polarization 
drift is 5 = v P i/cu = a/k v . For these reasons the mentioned 
gyro-center representation fails and the field magnitude is 
to be calculated at the actual position of the particle. An- 
other important feature is that v p i ~ k y , hence the stochas- 
tic heating is due to the electrostatic property of the wave. 
Although the resulting particle motion is deterministic, as 
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explained in iBellanl { 2006) the practical consequence of the 
described mechanism on the particle distribution function is 
the same as in an ordinary heating. 

The maximum achieved stochastic t emperature is 
l|McChesnev et al.lll987l ; ISanders et ai]|l998l ) 



T max = j(k 2 y p 2 t e<p/T l0 + 1.9) 2 ^. (6) 



The application of this effect to the heating of the solar 
corona by the ordinary density gradient dr iven drift wave has 
been pe rform ed in our recent pu b lication s IVranies fc Poedtsl 
(2009a) and IVranies fc Poedtsl <|2009bh . It is shown that 
the ions are more efficiently heated than electrons, and the 
heavier ions are heated better than light ions provided that 
kypi{e4>/Tio) 2 > 1.9. The nature of the heating is such that 
it acts mainly in the perpendicular direction, and it can also 
describe some other features of the coronal heating. 

The threshold potential (0 is presented in Fig. 6 in 
terms of X y , for the same parameters as earlier in the text. 
For the values above the curve, the stochastic heating takes 
place. On the other hand, as discussed earlier in the text, for 
large enough k z (f> the electric field exceeds the Dreicer value 
and the particle acceleration is in action too. The horizontal 
lines in Fig. 6 give k z (f> = Ed for the cases A z = 100 km and 
500 km. Hence, for the values of (j> above the both full and 
dotted lines, the plasma is subject to simultaneous heating 
of ions (in the perpendicular direction) , and an acceleration 
of bulk electrons in the parallel direction. 

The stochastic temperature (JSJ) is sensitive to the per- 
pendicular wavelength and the magnitude of the background 
magnetic field. Using the parameters from the previous text 
and the starting temperature of 1 MK one can easily ob- 
tain a temperature several orders of magnitude above the 
starting value, and this already for a very small amplitude 
of the potential <j>- The heating presented in Table 1 is for 
rather moderate values of the perturbed potential, i.e., for 
the left (lower) part of the curve in Fig. 6. From Table 1, 
one concludes that at very short perpendicular wave-lengths, 
the stochastic heating is in action already at very small am- 
plitudes of the wave potential. Also clear from Table 1 is 
the remarkable fact that for the same wave amplitude the 
heavier ions (helium) are heated more efficiently. 

The stochastic heating implies the condition (0 sat- 
isfied, regardless of the specific values of the two separate 
parts in that expression, while the instability analysis and 
the growth-time from Sec. 2 imply k y p 3 = k y pi < 1. For that 
reason we are formally allowed to estimate the growth time 
in Table 1 for X v ^ 2 m only, and this is of the order of a 
second, as shown in Sec. 3.1. For this wave-length the max- 
imum energy release rate per unit volume for the density 
no = 10 16 m -3 becomes 3noT max /(2fg) = 65 J/(m 3 s), and 
that is several orders of magnitude above the value accepted 
as necessary for a sustained heating in coronal loops and 
also well within the range of the total amount of the energy 
released in nano-flares. Some a dditional propertie s of the 
stoc hastic heating, presented in IVranies fc Poedtsl (|2009af ) 
and IVranies fe Poedtsl (|2009bh for the ordinary drift wave, 
remain valid for the present case as well. 



SUMMARY 

The heating of the solar corona by waves, and the propa- 
gation of waves in the solar environment has been the sub- 
ject of numerous studies in the past. Typically the wave 
heating model s rema i n within the widely used MHD, e.g. 
iPekiinlii et ail (|200ll ); ISuzukU (| 20041 ) . yet it is obvious that 
in such an app roach a lot of physics rem a ins out of scope as 
may b e seen in IVranies fc Poedtsl (|2006T ); |Pandev fc Wardlel 
(2008). The analysis performed in the present paper is aimed 
at showing even further some novel phenomena that follow 
from a multi-component description. It is based on the well 
established theory of the low frequency phenomena in mag- 
netized and inhomogeneous plasmas. It is crucially a multi- 
component plasma description. In the same time, it repre- 
sents a step forward in the rec ently published modeling o f 
the heating of the solar corona (|Vranies fc P oedts 2009allbl). 
The most important consequence of the here described fast 
growing mode is the electric field, with its drastically differ- 
ent scales in the parallel and perpendicular directions, that 
cannot be predicted within the widely used MHD theory. 
This electric field should be responsible for the acceleration 
of plasma particles and their simultaneous heating. The na- 
ture of the instability implies large scales in the direction 
of the magnetic field vector. In a realistic cylindric configu- 
ration the mode is weakly twisted around a magnetic loop, 
and such is the scene where the predicted acceleration and 
heating takes place. The essential difference between the in- 
stability presented here a nd the drift wave instability from 
IVranies fc Poedtsl (|2009al lbh is that the former can be well 
described within the multi-component fluid theory, while the 
latter is a purely kinetic effect. 

The gradient driven oscillations presented here are of 
relatively high frequency and therefore presently difficult to 
detect directly. Yet, they imply the presence of the electric 
fields that can be observed in coronal spectra due to line 
broad e ning a nd s hifts, that are measured and presented in 
iDavisI (|1977T ) and IZhang fc Smarttl (|l986T ). and attributed 
to plasma waves and instabilities, in particular to the lower- 
hybrid-drift and whistler instabilities. We show here that 
they can be described in terms of the gradient-driven drift 
instabilities. The possible consequences of the presence of 
gradient driven oscillations (the temperature anisotropy and 
a stronger heating of heaver ions) are well documented in 
the measurements of the spectral line widths of heavy ions 
ICranmer et all (|200SD . The mentioned strong E x B-plasma 
drifts should also be detectable as Doppler shifts in the spec- 
tra. 

One can with certainty claim that, at least in the start- 
ing stages of the presented gradient-driven instability, these 
processes are reasonably accurately described within the 
given model that can potentially be used for the descrip- 
tion of the particle acceleration and the heating of the solar 
corona. However, there exists a number of phenomena that 
will mostly negatively affect the proposed effects of heat- 
ing and acceleration, like the energy and particle diffusion 
and collisions, coupling to the Alfven wave, non-linearity, 
etc. Large wave amplitudes and the corresponding stronger 
heating will necessarily modify the starting plasma configu- 
ration, primarily the temperature gradient that is essential 
for the growth, and the density gradient too. Consequently, 
the values obtained here analytically may be far from ac- 
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curate. In particular, for large values of e<^/(«;Tj), various 
new effects must be included, and the model must be con- 
siderably improved in order to give accurate estimates in 
this limit. Though, the effective temperature presented in 
Table 1 goes over 300 million K, that is far above the tem- 
peratures that are needed. Hence, reducing the potential 
[i.e. the factor e(f> by two order of magnitude will still 
give desired temperature of a few million K. Such a reduc- 
tion of (j> will give lower values of the electric field k y <j>, yet 
this can be compensated by taking larger values of k v . For 
the given plasma parameters, the plasma Debye radius is 
around 0.5 mm, so that going to very short perpendicular 
wavelengths, even below the values from Table 1 is justified. 
Nevertheless, it is fair to say that a realistic (non-linear) de- 
velopment of the instability and of all the consequences that 
follow from it, can be described only numerically in codes 
that allow for a simultaneous change of the driving force 
(i.e., the inhomogeneous plasma background) in the process 
of the development of the instability. 
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